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We find three distinct phases; a tubular phase, a planar phase, and the spherical phase, in a 
triangulated fluid surface model. It is also found that these phases are separated by discontinu- 
ous transitions. The fluid surface model is investigated within the framework of the conventional 
curvature model by using the canonical Monte Carlo simulations with dynamical triangulations. 
The mechanical strength of the surface is given only by skeletons, and no two-dimensional bending 
energy is assumed in the Hamiltonian. The skeletons are composed of elastic linear-chains and rigid 
junctions and form a compartmentalized structure on the surface, and for this reason the vertices 
of triangles can diffuse freely only inside the compartments. As a consequence, an inhomogeneous 
structure is introduced in the model; the surface strength inside the compartments is different from 
the surface strength on the compartments. However, the rotational symmetry is not influenced by 
the elastic skeletons; there is no specific direction on the surface. In addition to the three phases 
mentioned above, a collapsed phase is expected to exist in the low bending rigidity regime that was 
not studied here. The inhomogeneous structure and the fluidity of vertices are considered to be the 
origin of such variety of phases. 

PACS numbers: 64.60.-i, 68.60.-p, 87.16.Dg 



I. INTRODUCTION 

A crumpling of surfaces has been investigated on the 
basis of the singularity analysis, and progress has been 
recently made on understanding the crumpling phenom- 
ena; the universal structure on the crumpled thin sheets 
was found in the formations of singularity of ridges and 
cones [H, . A similar transition to this phenomena was 
also found experimentally between the smooth state and 
the crumpled state in an artificial membrane, which is 
partly polymerized [3(. 

Studies have also been focused on the transition in the 
surface model of Helfrich, Polyakov and Klcinert (HPK) 
Jl,Ja, H from the viewpoint of statistical mechanics [7], 
SB EE El GJ, El- The bending rigidity is known to 
be stiffened by the thermal fluctuation of the surface, 
and this was confirmed in the statistical mechanics of 
membranes [3, [TBI, EH, E3, EH • Numerical studies were 
made to understand the transition in triangulated surface 

models [n H3, nn, m, m, m, m, m, H3, 121, m, n^i . The 

transition was reported as first-order in recent numerical 
studies HEHll. 

On the other hand, the concern with inho mog eneous 
surfaces has been growing over the past decade [33|,l3J|. A 
homogeneous artificial membrane that is coated by elas- 
tic skeletons is also considered to be an inhomogeneous 
membrane. Some of the mechanical properties of such 
membranes were revealed experimentally (3o| . The hop 
diffusion of membrane protein or lipids was observed, and 
as a consequence the compartment of cytoskeletons was 
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confirmed to be in biological membranes [36j |. It is also 
well known that the microtubule, which is an element 
of the cytoskeleton, gives a mechanical strength to the 
surface of the biological membranes. 

However, the surface collapsing phenomena and the 
surface fluctuation phenomena are almost unknown in 
such inhomogeneous models for membranes. Therefore, 
it is worthwhile to study an inhomogeneous fluid sur- 
face model within the framework of the conventional sur- 
face model of HPK. We note that the inhomogeneity in 
our model corresponds to the cytoskeletons in biologi- 
cal membranes as stated above. The fluidity realized by 
dynamical triangulations in the inhomogeneous model, 
as well as the fluidity in the homogeneous surface mod- 
els, corresponds to the lateral diffusion of lipids in mem- 
branes. 

In this paper we study a compartmentalized surface 
model by Monte Carlo (MC) simulations. The Hamil- 
tonian of the model includes no two-dimensional bend- 
ing energy but a one-dimensional bending energy. The 
model is defined on dynamically triangulated surfaces, 
where the free diffusion of vertices is confined inside the 
compartments. The mechanical strength of the surface is 
given only by the compartment boundary, which is com- 
posed of one-dimensional elastic chains and rigid junc- 
tions. Because the collapsed phase is expected to appear 
at sufficiently small bending rigidity b[kT] -^0(6^0), we 
concentrate on the phase structure at relatively large b 
in this paper. Consequently, information on the phase 
boundary at b— >0 remains unanswered. 

We recently reported numerical results of three types 
of surface models [371 . |38| | , which are similar to the model 
in this paper. Then, we should comment on the similar- 
ity/difference between the model in this paper and the 
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models in [37], [38j]. Firstly, the lattice structure of the 
model in this paper is very similar to that of the first 
model in [§3] and that of the model injpl], and is iden- 
tical to that of the second model in [37fl . Secondly, the 
lattice in this paper and that of the first model in [37| 
are the dynamically triangulated one, while the lattice of 
the second model in [33] and that in [38[ are the fixed- 
connectivity one. Thirdly, the Hamiltonian is different 
from the one in the first model in [37| . The Hamiltonian 
of the model of this paper includes only one-dimensional 
bending energy, which is defined on the compartment 
boundary while the Hamiltonian of the first in [37] in- 
cludes only a two-dimensional bending energy, which 
is defined all over the surface, and no one-dimensional 
bending energy is given to the compartment boundary. 
Therefore, the model in this paper is different from the 
three models in [13, HI . 

Our results obtained in this paper show that the model 
undergoes a first-order transition between the smooth 
phase and the crumpled phase. Moreover, the smooth 
phase can be divided into the spherical phase and the 
planar phase, and the crumpled phase can also be divided 
into the tubular phase and the collapsed phase, which is 
expected to appear at sufficiently small b because no self- 
avoiding property [H, EH is assumed in the model. 
It must be emphasized that such variety of phases can 
be seen neither in the conventional surface models nor in 
the compartmentalized models such as those in [37], IH[ • 

One remarkable result is the appearance of planar sur- 
faces. The echinocytic shapes of erythrocytes were ex- 
tensively studied, and they are currently known to be 
described by many models such as the area difference 
bilayer model [34j . The shape of membranes is also sen- 
sitive to the flow fields [H[ . Our model in this paper indi- 
cates that one possible origin of such planar shape comes 
from the inhomogeneity due to the cytoskeltal structure 
and the fluidity of lateral diffusion of vertices. 



II. MODEL 

Figure QJa) shows a triangulated surface of size 
(N,N s ,Nj,L) = (2322,600,42,6), where N is the total 
number of vertices including the junctions, N$ is the total 
number of vertices on the chains, Nj is the total number 
of junctions, and L is the length of chains between the 
two nearest-neighbor junctions. It should be noted again 
that Nj is included in N; junctions are counted in the 
total number of vertices. The junctions are assumed as 
rigid plates; twelve of them are pentagon and all the oth- 
ers are hexagon. The junction size in Figflja) is drawn 
many times larger than that of the lattices for the sim- 
ulations; and it will be discussed in the last part of this 
section. Thick lines on the surface in FigQJa) denote the 
chains, which are terminated at the junctions. 

The construction of the lattices is as follows: Let us 
start with the icosahedron. Every edge of the icosahe- 
dron is divided into £ pieces of uniform length, and then 




(a) (N, Ns,Nj, L) = (b) A rigid junction 

(2322, 600, 42, 6) with the chains 

FIG. 1: (Color online) (a) Starting configuration of surfaces 
of size (N, Ns, Nj,L) = (2322, 600, 42, 6), and (b) angles 
in the bending energy & of Eq.(f4]). Thick lines in (a) denote 
the compartment boundary composed of the linear chains and 
the rigid junctions of the hexagonal and the pentagonal plates, 
whose size is drawn many times larger than that of the lattices 
for the simulations. 

we have a triangulated surface of size No = 10£ 2 + 2 (= 
the total number of vertices on the surface). The com- 
partmentalized structures are constructed by dividing £ 
further into m pieces (m = 1, 2, • • • ). Thus, we have the 
chains of uniform length L — (£/m)—2 when m divides 
I. The reason for the subtraction —2 is because of the 
junctions at the two end points of the chain. Because 
the compartmentalized structure is a sublattice, the to- 
tal number of junctions Nj is given by Nj — 10m 2 + 2. 
The total number of bonds in the sublattice is 3Nj — 6, 
and each bond contains L — l vertices, then Ns is given 
by Ns = (3Nj — 6)(L — 1), which can be written as Ns — 
30m(£—3m). The hexagonal (pentagonal) rigid junctions 
are composed of 7 (6) vertices, then iVj — 12 hexagonal 
rigid junctions and 12 pentagonal rigid junctions reduce 
the total number of vertices Nq by (Nj— 12)x6 and 12x5. 
Therefore, we have ^ = ^0 — 6^,7+12, which can also be 
written as N — 10£ 2 —60m 2 +2. The thermodynamic limit 
of our model is defined by N — > oo, Ns — + oo, and Nj — > oo 
under the condition that L is finite. We have the thermo- 
dynamic limit at £— >oo and m— >oo. The lattice of size 
(N, N s , Nj, L) = (2322, 600, 42, 6) in FigOJa) is given by 
two independent integers (I, m) = (16, 2). 

The surfaces can be characterized by the length L. In 
this paper, we assume three values for L such that 

L = 6, L = 8, L = ll. (1) 

The value of L has a one to one correspondence with the 
total number of vertices n in a compartment; in fact, the 
values of L in Eq. (fBcorrespond to n = 2l, n = 36, and 
n = 66, respectively 13711 . We note that the effective phys- 
ical meaning of increasing (decreasing) L can be consid- 
ered as the increasing (decreasing) temperature. In fact, 
the surface fluctuation mainly comes from the thermal 
fluctuation of vertices inside the compartments. Because 
no bending energy is assumed inside the compartments, 
the fluctuation of vertices becomes large not only in the 
in-plane directions (free diffusion) but also in the direc- 
tion perpendicular to the surface. Thus, we consider that 
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the fluctuations are expected to grow with increasing n, 
i.e., increasing L. 

We use the surfaces of size (N,Ng,Nj) listed in Ta- 
ble |TJ Three different sizes (N, N$, Nj) are assumed for 
each L. The corresponding integers (I, m) are as follows: 
(16,2), (24,3), and (32,4) for the L = 6 surfaces, (10,1), 
(20,2), and (30,3) for the L = 8 surfaces, and (13,1), 
(26,2), and (39,3) for the L=ll surfaces. 

TABLE I: The surface size assumed in the simulations. Three 
sizes (N, Ns, Nj) are assumed for each L. 

L (N,N s ,Nj) (N,N s ,Nj) (N,N s ,Nj) 

6 (2322,600,42) (5222,1350,92) (9282,2400,162) 

8 (942,210,12) (3762,840,42) (8462,1890,92) 

11 (1632,300,12) (6522,1200,42) (14672,2700,92) 

The model is defined by the partition function 
,/ N 

Z = Y^ / l[dX i exp[-S(X,T)], (2) 
r ^ i=i 

S(X,T) = S 1 + bS 2 , 

where Si is the Gaussian bond potential, which is de- 
fined all over the surface, and S 2 is the one-dimensional 
bending energy, which is defined on the compartment 
boundary and will be given below. The parameter b is 
the bending rigidity. The integration symbol J in Eq. |(2J) 
denotes that the center of mass of the surface is fixed. 
YIt denotes the sum over all possible triangulations 7~, 
which are performed by the bond flip technique keeping 
the compartments unflipped. The bond flip procedure 
will be given in the following section. 

The integration measure ]X_i dXi ^ s gi ven by the 
product 

n In' \ / Nj \ 

j[dXi= [i[dx iq ?\ n^n^w ■ 

i=l \i=l J \i=l j(i) J 

(a = 3/2, 0), (3) 

where N' (=N — Nj) is the total number of vertices ex- 
cluding the junctions, Yli=i dXiqf denotes the integra- 
tion over the 3-D translational degrees of freedom (DOF) 
of the vertices i, and Yl%=i dXi Ylj(i) I'jU) denotes those 
of the 3D translational DOF and the 3D rotational DOF 
of the junctions i. The co-ordination number qi is the 
total number of bonds meeting at the vertex i, and qju\ 
is the total number of bonds meeting at the corner 
of the junction i. 

The parameter a was chosen to be a = 3/2 in [43l. |44|. 
while a = in many previous simulations on dynami- 
cally triangulated surfaces in the literatures. It is easy to 
understand that large positive a suppresses the configu- 
rations with large coordination number. Therefore, it is 
interesting to see the dependence of the phase structure 
on a. 



We chose both a = 3/2 and a = for the weight 
qf [H, HI] , and see whether the phase structure of the 
model depends on a or not. If the parameter is chosen 
to a = 3/2, then the coordination number qi serves as a 
weight of the integration dXi , while a = gives the uni- 
form weight. The weight YIa—i qf can also be written 

as niLi9? = exp(aX\ lo S?i)i and therefore, T^Ii 1? 
is considered to be the co-ordination dependent term 
~ a Ei'°S?! m the Hamiltonian; — a: log (ft changes 
its value only on dynamically triangulated surfaces. 

The Gaussian term Si and the bending energy term 
S2 are defined by 

Sx^iXi-Xj) 2 , S 2 = [1- cos 9 m ], (4) 
(H) (ij) 

where 53 (y) m ^1 ^ s the sum over bonds (ij) connecting 
the vertices i and j, and Y^Uj) m ^2 is also the sum over 
bonds (ij). in S2 is the angle between the bonds 

i and j, which include virtual bonds. The virtual bonds 
denote the lines between the center and the corners of the 
junction; the hexagonal (pentagonal) junction contains 
six (five) virtual bonds. 

Figure [TJb) is a junction and the chains linked to the 
junction on a fluctuating surface. Triangles are elimi- 
nated from the figure. One 0(ij) shown at a corner of 
the junction is defined by using a virtual bond and a real 
bond in a chain, and the other Quj) shown at a vertex is 
defined by real bonds on the same chain. 

The size of the junctions can be characterized by the 
edge length R, which is fixed to 

R = 0.1 (edge length of the junctions). (5) 

The value R = 0.1 is relatively smaller than the mean 
bond length 0.707, which corresponds to the relation 
Si/N — 1.5 satisfied in the equilibrium configuration of 
surfaces without the rigid junctions. As we will see later, 
the relation S±/N = 1.5 is slightly violated in the model 
of this paper because of the rigid junctions. 

Here we comment on the unit of physical quantities. 
Let a be the length scale of the model, then the unit 
of physical quantity that has the length unit can be ex- 
pressed by a; the unit of Si is [a 2 ]. The surface tension 
coefficient A in XSi+bS2 has the unit [kT/a 2 ] and as- 
sumed to be A= l[kT/a 2 }, and the bending rigidity b has 
the unit of [kT] as described above. 

Note that the bending rigidity b in the Hamiltonian is 
a microscopic quantity from the view point of statistical 
mechanical model, and therefore b is not always identical 
to the macroscopic bending rigidity of real physical mem- 
branes. However, the microscopic value b of real mem- 
branes can effectively be varied with the temperature, 
because b has the unit of kT. Therefore, it is possible to 
consider that the phase structure described in terms of b 
in the surface model corresponds to the phase structure 
described in terms of T in real physical membranes. The 
length scale a in the model is also a microscopic quantity 
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and, we consider that a is sufficiently smaller than the 
membrane size. 



III. MONTE CARLO TECHNIQUE 

A sequence of random numbers called Mersenne 
Twister [45| is used in the canonical MC simulations. 
The Metropolis technique is applied to update A and T, 
where the variable A denotes the position of the ver- 
tices and that of the junctions. The vertex position 
X is shifted so that X' = X + 5X, where 5 X is ran- 
domly chosen in a small sphere. The new position X' is 
accepted with the probability Min[l, exp(— AS*)], where 
AS = S'(new) — S'(old). The position A of a hexagonal 
(or pentagonal) junction, which is not a point but a rigid 
plate, is also integrated out by performing 3D random 
translations and 31? random rotations. 

Thus, the variable A is updated by a random N' (= 
N — Nj) shifts of vertices, a random Nj translations of 
junctions, and a random Nj rotations of junctions. These 
updates are denoted by (N' , Nj, Nj) updates of A. The 
N' shifts of A can be divided into N$ shift of the vertices 
on the linear chains and N' — N$ shifts of all the other 
vertices, which are those inside the compartments. 

The radius of the small sphere for 8X is fixed at the 
beginning of the MC simulations in order to maintain 
about 50% acceptance rate. The vertices on the linear 
chains carry the bending energy S2 in Eq.(|3]), while all the 
other vertices inside the compartments does not. There- 
fore, the acceptance rate is independently controlled in 
the two-groups of vertices. The radius for the random 
translation of the junctions and that for the random rota- 
tion are also independently chosen so that the acceptance 
rates are both about 50%. 

The summation over T in Z of Eq. (J2J) is performed by 
using the standard bond flip technique [22j, |23| . The flip 
is accepted with the probability Min[l, exp(— AS)]. The 
acceptance rate for the bond flip is not under control and 
is about 75%, which is almost independent of b. 

The bonds are labeled with sequential numbers. The 
total number of bonds is denoted by N' B , which excludes 
the number of bonds on the linear chains because the 
bonds on the linear chains remain unflipped. 

The bond flip is performed as follows: Firstly, the odd- 
numbered bonds are sequentially chosen to be flipped for 
the N' B /2 updates of T, and after that, the (A, Nj,Nj) 
updates of A are performed. Secondly, the remaining 
even-numbered bonds are chosen to be flipped for the 
N' B /2 updates of T, and after that, the (A, Nj, Nj) up- 
dates of A are performed. Thus, the (A, Nj, Nj) updates 
of A and the N' B /2 updates of T are consecutively per- 
formed, and these make one MCS (Monte Carlo Sweep). 

We introduce the lower bound 1 x 10~ 8 to the area of 
triangles. No lower bound is imposed on the bond length. 



IV. RESULTS OF SIMULATION 

A. a = 3/2 

As mentioned in Section [Til we assume the value of 
a in Eq.([3]) as a = 3/2 and a = 0. In this subsection, 
we present the results obtained under a — 3/2 by using 
snapshots and figures, and in the next subsection we will 
show some of the results under a = 0. 

The thermalization MCS is 1 x 10 7 in almost all cases. 
However, more than 1 x 10 s thermalization MCS were 
done close to the transition point in such cases that the 
surface is trapped in one phase at first and then changes 
its phase to a more stable one under a given condition. 
The total number of MCS for the production of samples 
is 0.8 x 10 8 ~ 1.3 x 10 s . At the transition point, about 
2 x 10 8 MCS was performed after the thermalization in 
some cases. 



(a) b = 21.2 (b) b= 21.8 (c) b= 22 




(d) The section (e) The section (f) The section 



FIG. 2: (Color online) The snapshots of surfaces of size 
(N,Ns,Nj,L) = (8462, 1890,92,8) obtained at (a) b = 21.2 
(tubular phase), (b) 6 = 21.8 (planar phase), and (c) b — 22 
(spherical phase), and (d),(e),(f) are the surface sections of 
(a),(b),(c), respectively. a = 3/2. 

We show snapshots of the (N,Ng,Nj,L) — 
(8462,1890,92,8) surface in Figs^a)^). They were 
obtained at (a) b = 21.2, (b) b = 21.8, and (c) b = 22, 
which respectively corresponds to the tubular phase, the 
planar phase, and the spherical phase. The snapshot of 
FigIU(b) at 6 = 21.8 was the final configuration produced 
after 2 x 10 8 MCS including 1 x 10 8 thermalizaion MCS; 
the planar surface was stable after the thermalization 
MCS. The surface sections are shown in FigsJ21(d)-[5{f); 
the sections in Figsl^d) and^e) were obtained by slicing 
the surfaces perpendicular to the vertical axis, and the 
section in FigslDJf) was obtained by slicing the surface 
perpendicular to a horizontal axis. All of the snapshots 
were drawn in the same scale. The axis of the tubular 
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surface Figj2][a) as well as the axis perpendicular to the 
planar surface FigJ^b) is spontaneously chosen. 

The planar phase is stable only on the L — 8 surfaces, 
while it seems unstable on the L = 6 surfaces and on the 
L = 11 surfaces. Even if the planar phase once appears on 
the surfaces of L = 6 and L=ll of size at least iV<9282 
and N < 14672, respectively, it eventually collapses into 
the tubular phase. Therefore, we find that no planar 
phase can be seen on the L — 6 and the L — ll surfaces; 
the tubular phase and the spherical phase are connected 
by a discontinuous transition on those surfaces. Thus, 
we understand that the planar phase appears depending 
on the size of the compartments. We should note that 
the planar surface may bend and fluctuate in the limit of 
N —>oo, and the tubular surface may also bend and wind 
in the same limit. 
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FIG. 3: The one-dimensional bending energy S2/N' s against 
6 obtained on the surfaces of (a) L = 6, (b) L = 8, and (c) 
L — 11. N' s (= Ns + 6Nj - 12) is the total number of vertices 
where S2 is defined. 

Figures EIa)|3]Jb), and (3Jc) show the bending energy 
Sa/N's °f Eq.(TJ]) against 6, which were obtained on the 
surfaces of L = 6, L = 8, and L = 11, respectively. 
A r g(= Ns+QNj— 12) is the total number of vertices where 
S2 is defined. 6Nj — 12 is the total number of corners 
of the junctions, which include 12-pentagons. The solid 
lines on the data were drawn to guide the eyes. Dashed 
lines drawn vertically denote the phase boundary be- 
tween the tubular and the spherical phases, the bound- 
ary between the tubular and the planar phases, and the 
boundary between the planar and the spherical phases. 
The discontinuous change of S2/-/V5 between the tubu- 
lar phase and the spherical (or the planar) phase is very 
clear in the figures and considered to be a sign of the 
first-order transition. 

In order to see the difference between S^/iVg in those 
three phases, we plot in Figs|3Ja)rjJb), andfjjc) the varia- 
tion of S 2 /N' s against MCS obtained at 6 = 21.2, 6 = 21.4, 
and 6 = 21.8 on the {N, N s , Nj, L) = (8462, 1890, 92, 8) 
surface. The thermalization MCS were not discarded; 
they were included only in those variations. S^/N'g at 
6 = 21.2 in FigfJJa) shows a jump from the spherical phase 
to the planar phase and a jump from the planar phase to 
the tubular phase; the corresponding MCS at the jumps 
were indicated with the dashed vertical lines. We also 
find in Figf5](b) a jump from the spherical phase to the 
planar phase. A jump is also seen in S^/iVg at 6 = 21.8 
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FIG. 4: The variation of Sz/N' s against MCS, which were 
obtained on the (N, N s ,Nj, L) = (8462, 1890, 92, 8) surface at 
(a) 6 = 21.2, (b) 6 = 21.4, and (c) 6 = 21.8. The dashed lines 
denote the MCS where the jumps occurred. The correspond- 
ing normalized histogram h(S2) obtained at (d) 6 = 21.2, (e) 
6 = 21.4, and (f) 6 = 21.8. The parameter a was fixed to 
a = 3/2. 



in Figf?]Jc) from the spherical phase to the planar phase. 

The value of 6 = 21.2 corresponds to the tubular phase, 
whereas 6 = 21.4 and 6 = 21.8 correspond to the planar 
phase, because the final states are considered to be stable 
states. The surfaces at 6 = 21.2 and 6 = 21.8 can be seen 
in the snapshots in Figs[3]Ja) and[^b). 

The distribution of S2/N' s are shown as the normal- 
ized histograms h(S2) in Figs|4jd)-[Hf), which respec- 
tively correspond to the variations in Figs|3Ja)-[lJc). We 
see that h(S2) in FigfJJd) has three peaks; two of them 
are almost overlapping and the other one is distinctly 
separated from the previous two. Those three peaks in 
h(S2) correspond to the spherical phase, planar phase, 
and the tubular phase. Two almost overlapping peaks 
can also be seen in h(S-2) in Figs|4je) and[4jf), and they 
are corresponding to the spherical phase and the planar 
phase. We remark that the surfaces hardly change not 
only from the tubular phase to the smooth (= spherical 
or planar) phase but also from the planar phase to the 
spherical phase on the L = 8 and L = 11 surfaces. For 
this reason, we find in Figs|3Ja)-[IIc) no jump-back from 
a higher 52 state (such as the tubular state) to a lower 
S2 state (such as the planar state). 

The two-dimensional bending energy is defined by 



1 



n i)- 



(6) 



where is the unit normal vector of the triangle i, and 
rii nj is defined on the common bond (ij) of the triangles 
i and j. S3 is not included in the Hamiltonian and is 
defined even on the edges of the rigid junctions. Figures 
IH a )"l3 c ) show S3 /Nb against 6 obtained on the surfaces 
of L = 6, L = 8, and £, = 11, where Nb is the total number 
of bonds including the edges of the junctions. The jump 
of S^/Nb in Fig^b) is clearly seen between the tubular 
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FIG. 5: The two-dimensional bending energy S 3 /N B against FIG. 7: The mean square size X 2 against b obtained on the 
b obtained on the surfaces of (a) L = 6, (b) L = 8, and (c) surfaces of (a) L = 6, (b) L = 8, and (c) h = \\. 
L — ll. Nb is the total number of bonds where S3 is defined. 



phase and the planar phase. On the contrary, S^/Nb 
in the planar phase in Fig[^b), as well as S2/N' s in the 
planar phase in Fig|3jb) , is not so clearly distinguishable 
from that in the spherical phase. 



S,/N 
1.518 

1.516 



S,/N 



S./N 



< 

tubular 


r» - 

spherical 




■ 


N=9282 

L=6 

0=1.5 








tubular 


spherical 










N=6522 






1=11 






a=1.5 











(a) 



9 9.5 
b 



20 21 22 23 



(b) 



(c) 



46 48 50 52 
b 



FIG. 6: The Gaussian bond potential Si/N against b obtained 
on the surfaces of (a) L — 6, (b) L — 8, and (c) L = ll. 

It is expected that the Gaussian bond potential Si/N 
is influenced by the phase transitions. The potential 
Si/N should be Si/N ~ 3/2, which is satisfied in the 
model without the rigid junctions because of the scale 
invariant property of the partition function in that case. 
However, the junction size R in Eq.([5]) is finite in the 
model of this paper, and therefore Si/N can slightly de- 
viate from 3/2. 

Figures[(JJa)-[5]Jc) show S%/N against b obtained on the 
surfaces of (a) L = 6, (b) L = 8, and (c) L = ll. Discontin- 
uous changes in Si /N shown in the figures are consistent 
with the discontinuous transitions of the model, although 
the changes are very small compared to the value of Si /N 
itself. We find also the expected deviation of Si/N from 
3/2 in the figures. 

Figures[7Ja)-[7Jc) show the mean square size X 2 , which 
is defined by 



X 2 



X) 



X 



to 



(7) 



where X is the center of mass of the surface. We see 
that the phase transition is not reflected in X 2 on the 
L = 6 surfaces in FigJTJa), and the transition is also not 
reflected in X 2 on the L = 8 surfaces in Fig[7]Jb) at the 
transition point between the planar phase and the spher- 
ical phase. To the contrary, X 2 discontinuously changes 



in FigtTf b) at the transition point between the tubular 
phase and the planar phase and also at the transition 
point in FigJTJc). All of these behaviors of X 2 at the tran- 
sition points are consistent with those of S2/N' s , S3/NB, 
and Si/N. 



B. a = 

In this subsection, we present some of the results ob- 
tained under a = 0. 
(a) b=20.9 (b)b=21.4 (c) b= 21.8 




(d) The section (e) The section (f) The section 

FIG. 8: (Color online) The snapshots of surfaces of size 
(N,N s ,Nj,L) = (8462, 1890,92,8) obtained at (a) b = 20.9 
(tubular phase), (b) b = 21.4 (planar phase), and (c) 6 = 21.8 
(spherical phase), and (d),(e),(f) are the surface sections of 
(a),(b),(c), respectively. a~0. 



Snapshots of surfaces of a = are shown in FigsJSJa), 
13b), &[c), which respectively correspond to the tubu- 
lar phase (b = 20.9), the planar phase (b = 21.4), and 
the spherical phase (b = 21.8). The surface size is 
(N, N s , Nj, L) = (8462, 1890, 92, 8), which is identical to 
that in FigH The snapshot in FigE^b) at b = 21.4 is 
the final configuration produced after 1.9 x 10 8 MCS in- 
cluding 1 x 10 7 thermalizaion MCS; the planar surface 
was stable throughout the simulation. Thus, we find 
that three distinct phases are seen also in the surfaces 
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of L = 8, and that the planar phase is unstable on the 
surfaces of L — 6 and L — ll under the condition a = 0. 
Therefore, we consider that the phase structure of the 
model is independent of whether a = 3/2 or a = 0. 
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total number of junctions increases (decreases) with de- 
creasing (increasing) L. Therefore, the strength to or- 
der the surface increases (decreases) with decreasing (in- 
creasing) L. Then, we expect that the surface is ordered 
(disordered) at sufficiently small (large) L at given in- 
termediate value of b. Moreover, it seems possible that 
two competitive forces to order/disorder the surface are 
balanced with each other at intermediate values of L and 
consequently, some new phase appears depending on b 
at those L. Note also that the possibility of the appear- 
ance of planar phase is not completely eliminated on the 
surfaces of L = 6 and L = ll of sufficiently large size. 



FIG. 9: The one-dimensional bending energy S^/iVg against 6 
obtained on the surface of (a) L = 6, (b) L = 8, and (c) L—ll. 
N' s (=Ns + 6N,j — 12) is the total number of vertices where S2 
is defined. 

The one-dimensional bending energy S2/N' s obtained 
under a = is shown in Figs|9][a)-[9]Jc). A discontinuous 
change can be seen in S2/N' s not only in FigOJb) at the 
phase boundary between the tubular phase and the pla- 
nar phase but also in Fig^c) at the phase boundary be- 
tween the tubular phase and the spherical phase. A jump 
of S2/N' s in FigJHKb) at the transition point between the 
planar phase and the spherical phase is very small, and 
hence is hardly seen just the same as in FigJ3^b) under 
a = 3/2 in the previous subsection. Thus, we find no 
difference between S^/N'g of a = and that of a = 3/2. 
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FIG. 10: The mean square size X 2 against b obtained on the 
surfaces of (a) L = 6, (b) L — 8, and (c) L — ll. 

The mean square size X 2 are shown in Figs llOf a)- 
[TOT c). A jump is also seen in X 2 on the L — 8 and 
L = ll surfaces in FigsfTUTb) and [TUT c), and it is hardly 
seen on the L = 6 surfaces of size up to (N, N$, Nj, L) — 
(9282,2400,162,6). These results are identical to those 
observed in Figsl7ja)-[7^c) under a = 3/2. 

Finally, we comment on the planar phase appeared 
only on the L — 8 surface. The thermal fluctuation of 
vertices inside the compartments disorders the surface 
against the bending energy of the compartment bound- 
ary. Therefore, the strength to disorder the surface in- 
creases (decreases) with increasing (decreasing) L if N 
remains fixed, as stated in Section [TTJ On the other 
hand, the mechanical strength of the surface increases 
(decreases) with decreasing (increasing) L, because the 



V. SUMMARY AND CONCLUSION 

We have shown that a dynamically triangulated spher- 
ical surface has three distinct phases; the tubular phase, 
the planar phase, and the spherical phase, and that they 
are separated by discontinuous transitions. The first- 
order nature was very clear from the discontinuity in the 
bending energies S2 and S3 not only at the transition 
point between the tubular phase and the planar phase 
but also at the transition point between the tubular phase 
and the spherical phase. We know that the model has 
the collapsed phase at sufficiently small b, since the self- 
avoiding property is not assumed at least. Therefore, we 
expect that the model has four different phases including 
the collapsed phase, although the order of the transition 
between the collapsed phase and the tubular phase is un- 
known. 

The mechanical strength of the surface is given only 
by elastic linear-chains with rigid junctions. The 
triangulated surfaces are characterized by the size 
(N, 2V5, Nj, L), where N is the total number of vertices 
including the junctions, N$ is the total number of vertices 
on the chains, Nj is the total number of junctions, and L 
is the length of chains between the two nearest-neighbor 
junctions on the starting configurations. These four pa- 
rameters are not totally independent, because these are 
given by two independent integers {£, m), where m di- 
vides I. In fact, iV=10£ 2 -60m 2 -(-2, N s = 30m(£-3m), 
Nj = 10m 2 + 2, and L = (i/m)-2. 

We assumed three different values for L such that L = 
6, L = 8, and L=ll in the simulations. The edge length 
R of the rigid junction was fixed to be R = 0.1. The 
parameter a, which represents a weight for the three- 
dimensional integrations of the partition function, was 
assumed as a = 3/2 and a = 0. 

It is remarkable that the model has the planar phase, 
which is stable only on the surfaces with a specific struc- 
ture. In fact, the planar phase can be seen on the surfaces 
of L = 8, and it is unstable on the L = 6 and L = ll sur- 
faces. The planar phase appears in a narrow region on 
the 6-axis between the tubular phase and the spherical 
phase, and it is distinguishable from the spherical phase 
because a small but finite discontinuity can be seen in 
the bending energies S^/Afg and Ss/Nb- The gap of the 
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bending energy S2 at the transition point is very small, 
i.e., S% in the planar phase is almost identical to that 
in the spherical phase; however, the double peak struc- 
ture was clearly seen in the histogram of S2, which is 
included in the Hamiltonian. From this, we confirmed 
that the transition between the planar phase and the 
spherical phase is of first order. Our model in this pa- 
per indicates that one possible origin of planar shape of 
spherical membranes comes from the inhomogeneity due 
to the cytoskeltal structure and the fluidity of lateral dif- 
fusion of vertices. 

We have confirmed that the results obtained at a = 3/2 
in Eq.Q remain unchanged when a = 0. The phase 



structure of the fluid surface model in this paper is in- 
dependent of the choice of a at least for a — 3/2 and 
a = 0. Large scale simulations should be performed. It 
remains to be studied how large (£, m) are sufficient for 
the thermodynamic limit of the model. 
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